Asymmetric quantum error correcting codes 
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The noise in physical qubits is fundamentally asymmetric: in most devices, phase errors are much 
more probable than bit flips. We propose a quantum error correcting code which takes advantage 
of this asymmetry and shows good performance at a relatively small cost in redundancy, requiring 
less than a doubling of the number of physical qubits for error correction. 
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Introduction. The quest of a quantum computer 
has stimulated a lot of interesting developments 
in recent years. However, despite a remarkable 
progress, all physical devices realized so far do not 
allow to build even a very small computer. One cru- 
cial aspect is the noise control. Quantum computing 
faces two antogonist constraints: one should be able 
to manipulate and address the results of a computa- 
tion, and at the same time one must keep the noise 
level low. While some hardware architecture may 
help to achieve this compromise, it is clear now that 
there will never exist a quantum computer without 
efficient quantum error correction (QEC). 

The basic principles of QEC have been written 
down in 0, 0, HI Q, and a number of QEC codes 
have been developed since then [5|,|6( . However, most 
of them require in practice a high level of redun- 
dancy (in coding language, a low rate): the number 
of physical qubits needed to effectively protect one 
logical qubit is large. Generally one expects that 
a higher rate might be achievable with good codes 
when the length of the information block is large: 
in this limit, the uncorrected errors correspond to a 
correlated flip of a large number of physical bits, the 
probability of which gets exponentially small. The 
classical coding theory shows the existence of codes 
that become ideal (saturate the Shannon limit) when 
the size of the block tends to infinity. Furthermore, 
recent progress on so-called low density parity check 
(LDPC) codes has revived older ideas by Gallager 
(3 that produce efficient algorithms for a fast de- 
codingof codes with performance close to Shannon's 
limit H, 0, El • The generalization of these classical 
schemes for quantum error correction is made very 
difficult by the requirement that a quantum scheme 
should correct two types of errors: bit flips as well 
as phase errors. Another important constraint is the 
difficulty to perform operations concurrently on the 
same bit: an efficient error correction scheme should 
involve a relatively small (o(iV)) read-out operations 
on each bit. Here we propose a new family of QEC 



codes which work at relatively low redundancy (typ- 
ically < 2 — 3 physical qubits for one logical qubit), 
can correct many mistakes and allow parallelization. 
So far the main attempt at finding such codes is the 
work hj. It uses so-called self-dual codes which are 
taylored to deal with a noise which is symmetric in 
all channels. We argue that in the physical devices 
conceived so far, the noise is typically asymmetric 
(a phase error is much more probable than a bit 
flip) , and one can exploit this asymmetry to develop 
more efficient QEC codes. The construction that 
we propose here makes use of two standard classical 
codes which are among the most efficient ones: it 
handles the relatively rare bit errors through a Bose 
Chaudhuri Hocquenghem (BCH) code [12j and the 
more frequent phase errors through a LDPC code, 
with performances close to those of the most power- 
ful random LDPC codes [lflj. 

Physical noise The level of the noise in a single 
physical bit is conveniently characterized by the re- 
laxation time, T±, and dephasing time T2, the two 
parameters that enter Bloch equation for a single bit 
(spin) dynamics. Because the relaxation always im- 
plies dephasing, the dephasing rate I/T2 has a con- 
tribution from the relaxation processes and a pure 
dephasing: 1/T 2 = 1/(2T X ) +T 4> . Generally, there 
are many ways to control the relaxation rate: first, 
the relaxation between two states with energy dif- 
ference AE requires a transfer of energy to the en- 
vironment, the amplitude of which becomes smaller 
as AE — > 0. Furthermore, in many physical imple- 
mentations these two states are separated by a large 
barrier that makes transitions between them rare. 
The situation is completely different with the de- 
phasing rate which is physically due to the fluc- 
tuations of AE with time. All low frequency pro- 
cesses contributing to the AE(t) dependence result 
in the decrease of (exp(— i f AE(t)dt) correlator, i.e. 
lead to the dephasing. In this respect, a particularly 
damaging effect comes from omnipresent 1/f noise. 
Thus, it is not surprising that in almost all devices 



studied so far, the relaxation rate can be made much 
slower than the dephasing: in a ty pical NMR device 
Ti ~ 10 - 100s while T 2 ~ Is HfJ, in superconduct- 
ing phase qubits Ti ~ Wfxs while T 2 ~ 100ns 
in superconducting charge qubit T\ ~ 100ns while 
T 2 ~ Ins 0| and finally for spin dots T\ ~ 1/is 
while T 2 ~ 10ns [3113. 

In the following we shall therefore assume that 
in physical qubits the noise is strongly asymmet- 
ric. Specifically, we study a noise channel defined 
as follows. Noise acts independently on each bit. 
It induces a bit flip with probability p x , and in- 
dependently it induces a phase flip with probabil- 
ity p z . The original state of the system, IV'o)) 
is thus changed to \ifj) = JL [K)™' IV> > 
with probability pf'^fl - p z ) N -^> n *p¥ l (1 - 
Px) N ~ ^* m % where m,,nt € {0,1}. The chan- 
nel acts on bit i by applying an operator [7, G 
{I, o x ,o\,o x o l z }. 

CSS codes. Our family of codes is of the CSS 
type B- It consists of two independent encod- 
ing/decoding devices dealing separately with bit and 
phase flips, for a string of N physical qubits. It uses 
M z 'z-checks' and M x 'x-checks'. The a-th z-check 
is defined by a set V(a) G {1, . . . , N} and by the op- 
erator C* — Yli£v(a) a i- Similarly, the a-th z-check 
is defined by a set W(a) G {1, . . . , N} and by the 

operator = n i6H /(a) 

By construction, the z-checks and x-checks all 
commute with each other, and the original state 
V> ) is an eigenstate of all the operators C*,C X , 
with eigenvalue 1. As Ui either commutes or an- 
ticommutes with these check operators, the noise- 
perturbed state is an eigenstate of the operators 
C* , C%, . The decoding operation consists of three 
steps: (i) measure the eigenvalues of the check oper- 
ators, ii) infer from these eigenvalues what was the 
corrupting operator, (iii) apply the correction oper- 
ator. In more detail: 

Step (i): The a-th z-syndrom is defined as the 
number u a G {0, 1} such that C^\ip) = (1 - 2u a )\ip). 
Similarly, the a-th x-syndrom is defined as the num- 
ber v a G {0, 1} such that C£|V>) = (1 - 2« a ) l^- 
Step (ii): From the z-syndroms {u a }, a G 



{1, . . . , M z }, we compute N numbers {7 
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} 



such that, for each a G {1, . . . , M z }: J2iev{a) ~ 
u a (mod 2), with the smallest possible number of 
rn's equal to 1. From the x-syndroms {u Q }, a G 
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{1, . . . , M^}, we compute numbers {ni, . 
such that, for each a G {1, . . . , Mx}: ^ 
v a (mod 2), with the smallest possible number of n's 
equal to 1. 

Step (iii): generate = 

JliLi ( (T x)" i (°~l) mi IV')' If t nc error correction is 



successfull, one should find \tp') = \tp ). ^ 
A CSS code is thus characterized by the sets V(a) 
and W(a) defining the checks. In building such 
a code, one must ensure that all check operators 
commute. This imposes that Va G {1, . . . , M z }, 
Va' G {1,...,M X }, the cardinal of \V(a) U W{a')\ be 
even. It is useful to define the parity check matrices 
of the two codes. The matrix H z is a M z x N ma- 
trix with entries in {0, 1}, defined by = 1 if and 
only if i G V(a). Similarly, H x is the M x x N matrix 
defined by = 1 if and only if i G W(a). The com- 
mutativity condition is satisfied when H z (H X ) T = 
(using Boolean algebra, i.e. mod(2) additions). The 
z-codewords are strings of N bits Xi G {0,1} such 
that, Va, J2i-Hai x i = 0(mod 2). Any x-check a 
defines a z-codeword through Xj = 1 if i G W(a), 
and Xi = otherwise. Similarly, z-checks define x- 
codewords. Most of the research on QEC so far has 
focused on the design of relatively small codes with 
good distance properties. If for instance all pairs of 
x-codewords are at a Hamming distance > 2d+l, the 
code will correct any set of < t flip errors. While this 
suggest to build codes which maximize the smallest 
distance between codewords, this strategy is not nec- 
essarily optimal when dealing with large blocklength 
(N ^> 1). Instead, what is practically required is 
that the probability of an error is small and it turns 
out that the best classical codes often have (very 
rare) pairs of codewords which are pretty close to 
each other [l(| • We shall use this approach to con- 
struct our x-checks. 

z-checks: BCH code Our z-code is an efficient clas- 
sical construction, a binary primitive BCH code (see 
ref. ^3 f° r an extended presentation). The code 
depends on two parameter m, t. The first one deter- 
mines the Galois field GF(2 m ) which is used, and 
the number t is equal to the number of errors (phase 
flips) that the code can correct. The number of vari- 
ables (and therefore the number of qubits) is given 
by N = 2 rn — 1. If a is a primitive element of the 
field GF(2 m ), the powers a r , r G {1, . . . , N} are N 
distinct elements of the field, building a cyclic group 
under multiplication. At the same time, GF{2 m ) is 
a vector space of dimension m over GF{2): every el- 
ement a r can be decomposed as a r — J2™=o Irp ^ \ 
where the coefficients 7 are in {0, 1}. The check ma- 
trix H of the code is defined as 
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This matrix can be seen either as a t x N matrix 
with elements in GF(2 m ), but another interpreta- 



tion is also useful. If we write each element a r of 

H as the m component vector ... , we ob- 

V 7r(m-l) / 

tain the tm x TV parity check matrix H z with entries 
in GF(2) = {0, 1}. Therefore M z = tm. BCH de- 
coding relies on algebraic properties which are most 
easily written in terms of polynomials. Here we shall 
just present the basic result in the case t — 2. If two 
of the N bits are flipped by noise, and these indices 
correspond to the elements of GF(2 m ) called X , /3 2 , 
the check matrix H , applied to the error vector, 
gives two syndromes £i — /3 1 +0 2 an d C3 = 0\ + 02- 
Decoding consists in finding 11 2 given Ci;C2- It 
is easily seen that this system has a unique solution 
in GF(2 m ) (up to the permutation of 1 and 2 )- 
the code with t = 2 corrects exactly any set of < 2 
errors. The same construction works for arbitrary t, 
and good decoding algorithms exist: the code cor- 
rects any set of < t errors. In practice we have used 
the Berlekamp algorithm [l^, adapting some soft- 
ware available from (lif. 

Generation of the x checks: LDPC code. Some 
BCH codes are self-dual; in such a case one gets a 
quantum code using H x = H z [3]. But in order to 
get a much better performance (for large N) on the 
x-channel, we prefer to use a code as close as possi- 
ble to the random LDPC codes. The commutation of 
the x and z checks is obtained by the following pro- 
cedure. Given a BCH code with parameters m, t, we 
can generate a x-check a with any degree n > It + 1 
using a variant of its standard decoding algorithm. 
The first n — t elements of W(a) are chosen as a 
random subset of {1, ... , N} with distinct elements, 
taken uniformly among all such subsets. Let us call 
0!,..., n ^ t the corresponding elements of GF{2 m ). 
We look for the remaining t elements which are so- 
lutions of the decoding equations 



Vse{l,...,t}: J2(Pn- 



t+r) 



Provided that the solution of these equations ex- 
ists (which happens with probability 1/tl) the ele- 
ments I3 d _ t+1 . . . d can be found using any stan- 
dard BCH decoding algorithm, like Berlekamp's 
one. The indices corresponding to the elements 
0!, ... , n _ t ■ ■ ■ n form the subset W(a) defining 
the a-th x-check. As 1 , . . . n is a codeword of the 
BCH code, the commutativity condition is satisfied. 

Clearly, the indices in V{a) do not form a random 
subset of size n. However, if the map used in gen- 



erating I3 n _ t+1 ...0 n from lt 



is chaotic 



enough (we shall refer to this hypothesis as the 
'chaos hypothesis' in the following), one can hope to 
generate a set of x-checks with performances close 



to the ones of classical random LDPC codes. Tin's 
is what we have found numerically. In practice, for 
a given value of t, we generate a large enough pool 
of parity checks, all having degree n = It + 1. From 
this pool, we select a number M z of checks, in such 
a way that the degrees of the variables in the cor- 
responding factor graph has a narrow distribution. 
This is done by the following inductive procedure. 
At each step we order the remaining (unused) set 
of checks by their 'quality' which is defined as the 
number of minimal degree variables that would be 
affected by addition of this check. We then add one 
(randomly chosen) check of the highest 'quality' and 
repeat the procedure. 

The practical decoding of our LDPC code uses the 
standard 'belief propagation' (BP) algorithm|||j|, a 
message passing algorithm which is equivalent to an 
iterative solution of Bethe equations. 

Performance. An important parameter of the 
code is its degree of redundancy. We have checked 
that the various checks are generically linearly inde- 
pendent, so the z-rate (resp. x rate) is obtained as 
R z = 1 — (resp R x — 1 ~ -jf-) and the quantum 
rate of the code is R = 1 - M *+ M - . 

The error correction ability depends on the chan- 
nel. In the z-channel (bit flip errors), by construc- 
tion, the BCH code is able to decode up to t errors. 
Therefore the probability of error in decoding this 
channel is 
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pi z (l-p z ) N -i 



(3) 



which is well approximated, for the small values of p z 
which interest us here, by 1 — e~ Np * X]*-=o(-^P z ) : ' /j !- 
Let us now turn to the x-channel. The perfor- 
mance of BP decoding for random LDPC codes can 
■(2j)>e studied analytically in the limit of large block- 
length 10]. Within the chaos hypothesis, one could 
thus derive the threshold for zero error decoding in 
the large N limit. However in practice we are in- 
terested in not-too-large values of N. We have thus 
tested numerically the BP decoding of our x-code. 

The simulation is run as follows. We fix an 'ac- 
ceptable' value of the block error Pbiock for decoding 
N bits, both in the x and in the z-channel, in prac- 
tice Pbiock = 10~ 4 . For given values of N (or m) and 
t, eq.© gives the noise level p z that can be corrected 
in the z-channel, and the channel asymmetry gives 
the ratio p z /p x - We then test various x-codes, vary- 
ing M x until the block error in the x-channel is less 
than Pbiock- Results are summarized in the following 
table, which studies asymmetries p z /p x = 0.01, 0.1. 



Q_ 




sic low density of the code. We believe that these 
codes might be quite useful for the realistic physical 
implementation of quantum memory. We have not 
investigated the possibility of using them for fault 
tolerant quantum computation, this is the subject 
of the future research. 
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FIG. 1: +: block error in the :r-channel, Pg IT , versus the 
phase error probability p x , for the code with m = 12, 
t = 4, M x = 1378. The line is a guide to the eye. Also 
shown is the same curve for a random LDPC code (x). 
The inset gives the same data in a log-log plot 
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We see that large enough codes provide a good 
performance. For instance, m = 12, t = 6, code with 
N = 4095 qubits is able to correct a noise level of 
p z = 2.5 10~ 4 in the z channel and p x = 2.5 10~ 2 in 
the x channel with block error probability smaller 
than 10~ 4 . Notice that for these values of p z , p x , 
the probability of a block error without any error 
correction would be 1 — (1 — p z ,x) N , giving .63 for 
the z-channel and 1 for the x-channel. Figure 
gives the block error in the x-channel, Pg TT , versus 
the phase error probability p x , for one given code. 

Conclusions. We have provided an explicit con- 
struction of quantum codes with rates Q ~ 0.5 that 
are able to correct a few errors in one channel (bit 
flips) and have close to optimal peformance in an- 
other (phase errors), together with efficient decoding 
procedures. One important aspect of these codes is 
the fact that the number of operations to be done 
on one given bit is much smaller than N. In the 
z-channel this is due to the fact that we use a small 
value of t, in the x-channel it is due to the intrin- 
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